Asymmetric magnetic interference patterns in 0-n Josephson junctions. 
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We examine the magnetic interference patterns of Josephson junctions with a region of 0- and of 
7r-phase shift. Such junctions have recently been realized as c-axis YBCO-Pb junctions with a single 
twin boundary in YBCO. We show that in general the junction generates self-fields which introduces 
an asymmetry in the critical current under reversal of the magnetic field. Numerical calculations of 
these asymmetries indicate that they account well for unexplained features observed in single twin 
boundary junctions. 
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Josephson tunneling experiments on the high temper- 
ature superconducting compound YBa 2 Cu3C>7 (YBCO) 
have played an important role in establishing the pre- 
dominantly d— wave symmetry of the superconducting or- 
der parameter However, YBCO has orthorhombic 
symmetry which implies that the order parameter cannot 
be purely d-w&ve and must contain an additional s-wave 
contribution [p}-p~3|| . This has significant consequences for 
c-axis Josephson tunneling experiments between YBCO 
and Pb (a standard s-wave superconductor). In such 
junctions the d-wave component is forbidden by symme- 
try to contribute to the lowest order Josephson current 
so that the observed current j|] is due solely to the s- 
wave component of the YBCO. There is a clear difference 
between untwinncd and twinned YBCO samples in the 
observed current. The current is considerably suppressed 
for the latter compared to that of the former. This can 
be understood if we assume that the d-w&ve component is 
essentially uniform while the s-wave component changes 
sign at each twin boundary, i.e. a twinned sample yields 
a c-axis junction with alternating 0- and 7r-phase shift 
regions. This alternation of sign (or phase) gives rise 
to destructive interference effects for the total Josephson 
current in heavily twinned samples fl(i|j . 

To gain more insight into the c-axis Josephson tun- 
neling results Kouznetsov et. al. (m) have built c-axis 
YBCO - Pb junctions that contain a single twin bound- 
ary and have measured its magnetic interference pattern. 
These elegant experiments show that the critical current 
as a function of the applied magnetic field, Z c (H e ), be- 
haves qualitatively like a junction separated into two re- 
gions, one with 0- and the other with 7r-phase shift (such 
a model naturally arises due to ir phase change of the 
s-wave component as the twin boundary is crossed). A 
noteworthy unexplained feature of the data in Ref. H] is 
an asymmetry between J c (H e ) and J c (— H e ) for fields in 
the junction plane that are neither along the twin bound- 
ary nor orthogonal to it. In the short junction limit (for 
which the junction dimensions are much smaller than the 



Josephson penetration depth) such an asymmetry implies 
that time reversal symmetry is broken in the junction 
when there are no applied currents or magnetic fields 
|l5f . For example such an asymmetry will occur if there 
is local broken time reversal symmetry due to a phase 
difference of it/2 between the s and d components at the 
twin boundary (far from the twin boundary this phase 
difference must be or tt) ]T^ , p^ |. 
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FIG. 1. Junction geometries considered in this article. The 
left figure shows the junction plane while the right figures 
depict the two geometries considered in this article. 

In this article we present another origin for such an 
asymmetry in junctions that are time reversal invari- 
ant. In particular, such asymmetries arise whenever par- 
ity symmetry is broken in these junctions. This parity 
breaking can occur in two ways: either the junction is 
not invariant under a parity transformation about the 
center of the junction (intrinsic parity breaking) or the 
current running from the leads into the junction does 
not obey I\qs — I\-ds where OS is the boundary of 
the junction (e.g the current does not enter and leave 
a rectangular junction through opposite edges) (extrin- 
sic parity breaking). In the context of the — 7r junction 
model for a single twin boundary the junction has in- 
trinsic parity symmetry whenever the twin boundary lies 
in the center of the junction (called a symmetric junc- 
tion after Kirtley et. al. |llf|). A numerical study for 
a symmetric — 7r junction with extrinsic broken parity 
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symmetry is shown to account well for the experimen- 
tally observed magnetic diffraction patterns. Numerical 
studies of asymmetric — ir junctions with extrinsic par- 
ity conserving boundary conditions also reveals asymme- 
tries that should be observable in experiments like those 
of Ref. (l4| . This discussion is also of relevance to corner 
junction flux modulation experiments ]5|,p|,p"7| . 

We consider here junctions with the two geometries 
shown in Fig.l. The geometry G2 breaks extrinsic parity 
symmetry while Gl does not. The coordinates perpen- 
dicular and parallel to the twin boundary are called x 
and y, respectively. The order parameter of the Pb has 
standard s-wave symmetry, ^pb = \^>pb\e l(p . For YBCO 
there are two components (* d ,* s ) = e iv '(\^ d \, e ia \^ s \) 
where Vf^ and ^ s correspond to the d- and s-wave com- 
ponents, respectively. The relative phase a is fixed by 
the bulk free energy of YBCO and is for one type of 
twin domain and ir for the other |Io|Jll| l. The intrinsic 
phase shift of the junction is determined by a. Note that 
the choice of a equal to or n is a matter of convention 
(freedom of gauge in the two superconductors). However, 
once it is fixed somewhere in the junction it is determined 
everywhere. 

The Josephson phase, the local phase difference 4> = 
ip — ip' over the junction follows the relation H(x) x z — 
(<I>o/27rcO(V</>) while H satisfies the Maxwell equation 
V x H = 47rj/c ($ = hc/2e) (18). For the lowest or- 
der Josephson coupling (between \& a and ^fpb) we can 
express j z as j c (x) sin[0 + a(x)] (j c (x) > 0). By restrict- 
ing a(x) to take on the values or n and allowing j c (x) 
to take on any value this relation for j z describes any 
time reversal invariant junction. We therefore use this 
form for more general considerations and then specify 
j c (x) and a(x) to describe the — tt model for the single 
twin boundary junction for numerical results. Together 
these relations lead to the Sine-Gordon equation describ- 
ing the spatial dependence of <f> throughout the bulk of 
the junction 



V 2 = A J 2 (x) sin [<?!> + a(> 



(1) 



where A,/(x) = [<&q /2ndj c (x)] 1 / 2 with the boundary con- 
ditions 



Ajn • V0|as = n • (h x z)\ dS 



(2) 



where n is the outward normal to the boundary dS. 
The variable h is a function of the reduced variables 
h e _ 2^Ai He ( He is the a ppii ec j external field) and 

where j c is the average value of j c (x) and 
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A j is determined from j c . For the geometries shown in 
Fig 1 Ampere's Law implies that we have for Gl 
(extrinsic parity conserving boundary conditions) 



n • [h x z]\as = 



h e y +i dS = Bl 

-h% dS = B2 

-h% +i dS = B3 

hi dS = B4. 



(3) 



and for G2 (extrinsic parity violating boundary condi- 
tions) 



fh x z] 



dS 



h e y dS = Bl 
-h% +i dS = B2 
+i dS = B3 



(4) 



i><: 



as = BA. 



The Sine-Gordon equation and the boundary conditions 
can be recast into a variational form with the functional 



F = \jj s d 2 x (Vcj)) 2 j2 - A7 2 (x) cos^ + a(x)) 
+ f gs dx(f)[n- (h x z)]. 



(5) 



The positive and negative critical currents of the junc- 
tion are defined as the maximal positive and negative 
supercurrents, respectively, that can pass through the 
junction. They satisfy the following symmetry relations 
based on the structure of F: 

7 c [H e , a(x), Aj(x)] = - J c [-H e , -a(x), Aj(x)] (6) 
J c [H e , a(x), A, 7 (x)] = J c [H e , a(x) + a , Aj(x)] (7) 

where ao is a constant phase shift. Eqs. ^| and |^ follow 
from the time reversal and gauge invariance of the junc- 
tion respectively. For geometry Gl, due to the extrinsic 
parity symmetry, the following relation also holds 

7 c [H e , o(x), A. 7 (x)] = 7 C [-H e , a(-x), Aj(-x)] (8) 

more specifically this is a result of the invariance of F 
in Eq. (^) under the transformation h e — > — h e , a(x) — > 
a(— x), A(x) — > A(— x), and </>(x) — > <p(— x) (a product of 
time reversal and parity symmetries). No such relation 
holds for geometry G2 which has broken extrinsic parity 
(this is true for any geometry that exhibits extrinsic par- 
ity breaking). If the junction is time reversal invariant 
then the phase a(x) is restricted to be or 7r. For such a 
junction 7 c (H e ) = J c (— H e ) is required if the boundaries 
have extrinsic parity symmetry, if j c (x) = j c {— x), and if 
a(x) = a(— x) or if a(x) = a(— x) + ir. Consequently the 
observation of 7 c (H e ) ^ I c (— H e ) in junctions that are 
time reversal invariant is due to broken parity symme- 
try in the junction. In the following we will consider the 
— 7r model for the single twin boundary junction where 
Aj(x) = Aj and a(x) = irO(x — M) where O is the step 
function and —L/2 < M < +L/2 (M is the position of 
the twin boundary according to Fig.l). 

Consider the geometry Gl. The one dimensional model 
that results when the magnetic field is applied along the 
twin boundary has been previously studied by several 
groups (e.g. see Ref. (lqjl7 |). Using the relations (0) and 
(||) with ao = tt it is easy to see that J c (H e ) = 7 C (— H e ) 
if the twin boundary lies in the center of the junction 
(a symmetric junction |]l6|| ). For other positions of the 
twin boundary there is no symmetry constraint to enforce 
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2c (H e ) ^ ic(— H e ), however it can be shown for these 
asymmetric — tt junctions that 

h{h%,hl)=h{-h%,hl). (9) 

Numerical analysis shows 7 c (H e ) ^ J c (— H e ) when ft,^ 7^ 
for asymmetric — 7r junctions (see below). As pointed 
out earlier p7| the change of a from to 7r introduces a 
spontaneous flux line corresponding to a u tt vortex" , i.e. 
vortex of half the conventional phase winding and a flux 
$o/2, if the junction is much longer than the screening 
length Xj Q ■ Even for junctions with a length compara- 
ble to Aj such self fields appear. However, in this case the 
screening is imperfect and new effects due to the position 
of the twin boundary arise. For an asymmetric junction 
the screening is more effective for the longer of the a = 
or a = tt regions leading to J c (H e ) ^ J c (— H e ). Note that 
the time reversal invariance of asymmetric — tt junc- 
tions requires 7 c (H e ) = — I c {— H e ) so that experiments 
in which the positive and negative critical currents are 
averaged (such as in Ref. ||]) will not reveal this asym- 
metry. A recent careful numerical study of this model 
when the magnetic field is along the twin boundary by 
Kirtley et. al. |L6| only considered one sign of the applied 
field and consequently did not uncover this asymmetry. 
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FIG. 2. I c (H e )/j c L 2 for geometry G1,M/L = 0.25, and 
Xj/L =0.3, 0.5, and 1.0 from top to bottom respectively. 
Ic(He) represents the positive critical current and the negative 
critical current is given by —I c (—H e ). 

To examine the magnitude of the difference between 
lc(H e ) and J c (— H e ) for asymmetric — 7r junctions a 
numerical calculation is required. Since the two dimen- 
sional problem is numerically intensive a detailed study 
for the field along the twin boundary (for which the 
problem becomes one dimensional) was done. The tech- 
nique we use is to discretize the variable <f>(x) into N 
variables and to minimize the free energy by a quasi- 
Newton technique. The critical current was found by 
increasing the current until no stable minimum can be 
found anymore. The criterion for failure to reach a 
stable state in F in Eq. a is similar to that used by 



Kirtley et. al. (T(| and is given by e > 10 4 where 
(c/tt) 2 = [fa - 0a - (t - h)/a x / 2 ] 2 + [0at - 4> N ^ - (» + 
h)/a l ' 2 f + Et^+i + 4>i-i ~ Hi - sin(& + 9i)/a} 2 
where a — NX 2 /L 2 . Our numerical results using N — 
100 for the symmetric 0— 7r junction and for the 0—0 junc- 
tion agree well with those found by Kirtley et. al. [ fl6| . 
Fig. 2 shows the results for M / L = 1 /4 and for a variety 
of penetration depths. The asymmetry between I c (H e ) 
and I c (—H e ) is clearly visible and this further implies 
that zero field does not correspond to either a maximum 
or a minimum of the magnetic diffraction pattern. 

Now consider the geometry G2. This geometry is the 
relevant one for the experiments on single twin bound- 
ary junctions ]l4| ] and in this case there are no symme- 
try relations enforcing 7 c (H e ) = I c (— H e ) for symmet- 
ric — tt junctions due to extrinsic parity breaking. To 
gain an understanding of the effects of this geometry on 
the I c (H e ) pattern it is helpful to first consider — 
junctions (i.e. Aj(x) = Xj and a(x) = independent 
of x). Symmetry relations enforce 7 c (H e ) = I c {— H e ) 
when the field is along the dotted diagonal shown in Fig. 
1. These relations also imply I c (H e ,e) = I c (—H e ,—e) 
where e is defined in Fig. 1. A 2D generalization (for 
a 30 by 30 system and L/Xj = 5) of the ID numeri- 
cal method presented earlier shows that as expected a 
conventional 7 c (H e ) for e = arises and as e increases 
the central peak tilts so that the maximum I c moves 
from H e = to positive fields. The maximum I c is 
furthest away from H e = for e = tt/2. We have not 
been able to find any experimental reports of this angu- 
lar dependent J c (H e ) pattern for conventional Josephson 
junctions. To some degree this angular asymmetry still 
occurs for the — tt junction. For the — tt junction 
the interplay of the presence of the — it phase shift and 
the geometry makes a qualitative discussion difficult. In 
Fig. 3 we present numerical results on a 30 by 30 system 
for A j = 1/2. The numerical results exhibit all of the 
main features observed in experiments: the asymmetry 
is not apparent for the field along or perpendicular to the 
twin boundary but is so for other field orientations, the 
sign of 7 c (H e ) — I c {— H e ) changes when the field passes 
through the twin boundary, and the magnitude of the 
asymmetry is in good agreement with the experimen- 
tally quoted Xj/L values. Note that the asymmetry is 
large for the 8 = ir/4 and small for 9 = — it /A (in fact 
it appears to vanish between 9 = — tt /A and 8 = — tt/3) 
as was the case for the — junction. This agreement 
between experiment and theory gives further support to 
the claim that the s-wave component in the YBCO super- 
conductor changes sign as the twin boundary is crossed. 
These results are for a 30 by 30 system and studies for 
the ID system show that increasing the system size does 
not change the qualitative behavior of the J c (H e ) pat- 
terns but does increase the size of the critical currents 
and also moves the minima in the J c (H e ) patterns to- 
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wards zero field. This explanation for the experimental 
results has the consequence that 7 c (H e ) = — I c {— H e ) as 
follows from the time reversal symmetry of the junction. 
A deviation from this equality implies that time reversal 
symmetry is broken in the junction. It is also of inter- 
est to study experimentally the geometry Gl since some 
qualitative differences in the I c (H e ) patterns between the 
two geometries are expected to arise. For example in ge- 
ometry Gl the asymmetry should vanish for symmetric 
— 7r junctions and for asymmetric junctions as h% is re- 
versed 7 c (H e ) — I c (— H e ) should be the same due to Eq. || 
(in contrast to changing sign as it does for geometry G2). 
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FIG. 3. 7 c (H e ) I j c L 2 for geometry G2, a symmetric — tv 
junction, and Xj/L = 0.5. From top to bottom correspond 
to 8 — 0, ±7r/6, ±7r/4, ±7r/3, 7r/2 where positive (negative) 9 
is on the right (left) side. Successive plots have been offset by 
0.16. 

In conclusion we have shown that parity breaking in 
time reversal invariant Josephson junctions lead to an 
asymmetry in the critical current as the magnetic field is 
reversed. A detailed examination of the magnetic diffrac- 
tion patterns of — it junctions has shown that self-fields 
account well for unexplained features observed in single 
twin boundary junctions. Experimental tests have been 
proposed to further examine this theory. 
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